Quantum toy model for black-hole back-reaction 
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We propose a simple quantum field theoretical toy model for black hole evaporation and study 
the back-reaction of Hawking radiation onto the classical background. It turns out that the horizon 
is also "pushed back" in this situation (i.e., the interior region shrinks) but this back-reaction is 
not caused by energy conservation but by momentum balance. The effective heat capacity and the 
induced entropy variation can have both signs - depending on the parameters of the model. 
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Introduction Black holes are arguably the most sim- 
ple and at the same time most intriguing objects in the 
universe. The no-hair theorem states that they can fully 
be described by a small set of parameters such as their 
mass M and angular momentum J. Yet our standard 
picture of black holes contains many striking properties: 
Even though black holes should be completely black clas- 
sically, they emit Hawking radiation due to quantum ef- 
fects This evaporation process causes the black hole 
(horizon) to shrink (in the absence of infalling matter) 
due to the back-reaction of Hawking radiation. There- 
fore, black holes possess a negative heat capacity Q, i.e., 
the temperature grows with decreasing energy. Extrapo- 
lating this picture till the final stages of the evaporation, 
the black hole should end up in an explosion, where its 
temperature blows up and thus effects of quantum grav- 
ity should become important. Perhaps most fascinating 
is the observation that the second law of thermodynam- 
ics apparently Q requires to assign an entropy S to the 
black hole, which is determined by the horizon surface 
area A via S = A/ A (in natural units h= G = c = 1). 

Taking the analogy between black holes and thermo- 
dynamics seriously provides a very consistent picture, 
which has been confirmed by various gedanken experi- 
ments [3, 0] considering the construction of heat engines 
with black holes etc. It almost seems as if nature was try- 
ing to give us some hints regarding the underlying struc- 
ture which unifies quantum theory and gravity - which 
we do not fully understand yet. In order to understand 
these hints better, it might be useful to ask the question 
of whether (and how) the aforementioned properties de- 
pend on the detailed structure of the Einstein equations 
or whether they are more universal. For example, the 
study of condensed-matter based black hole analogues 
|, 1 shows that Hawking radiation is a fairly robust 
quantum phenomenom which just requires the oc- 
currence of an effective horizon and is quite independent 
of the Einstein equations. In contrast, the introduction 
of a black hole entropy with the desired properties seems 
to rely on the Einstein equations. 

In the following, we try to further disentangle universal 
features from properties which are specific to black holes 



(e.g., Einstein equations, rotational symmetry, conserved 
ADM mass). To this end, we propose a toy model which 
captures some of the relevant features of black holes and 
allows us to study the back-reaction of the emitted Hawk- 
ing radiation onto the classical background solution. 

Toy Model In the toy model we are going to discuss, 
the gravitational field will be represented by a real scalar 
field ip in 1+1 dimensions with the Lagrangian (h = 1) 



ci[d x ^Y)-v^) 



(i) 



With respect to the propagation speed c$ of the ip field, 
this form is Lorentz invariant. The potential V{tji) is 
supposed to be very stiff, i.e., the field ip is assumed to 
be heavy in the sense that it can be well approximated 
by a classical field. For definiteness, we choose the sine- 
Gordon potential V{ip) cx 1 — cos{ip / 'ipo) , but other po- 
tentials admitting stable solitonic solutions would also 
work. The global ground state "0 = then corresponds 
to a vanishing gravitational field whereas a kink (topo- 
logical defect) models a black (or white) hole horizon 

tp(x) = -40o arctan (exp {-£[x - x k ink]}) ■ (2) 

The position x — Xkink of the kink at rest is arbitrary and 
its width l/£ is determined by V(ip) and c^. In compar- 
ison to other models of black holes (see, e.g., [!,[§]), the 
advantage of the above set-up lies in the topologically 
protected stability and localization of the kink, which 
behaves very similar to a particle (see also floj|). 

In order to study Hawking radiation and its impact on 
the kink, we consider a massless quantum field <p coupled 
to the heavy field ip via the coupling constant g 



C^ = \([d t cp + giPd x cP] 2 -cl[d x cPf 



(3) 



Note that the velocity of the light (massless) field may 
differ from c^. The propagation of the light field <p in 
the approximately classical background ip is completely 
analogous to that in a gravitational field described by the 
Painleve-Gullstrand-Lemaitre metric (cf. [5, 6|) 



ds l 



[c\ - v 2 ) dt 2 - 2v dt dx - dx 2 



(4) 
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where v = gip denotes the local velocity of freely falling 
frames. A horizon occurs if this velocity v exceeds the 
speed of light c^. Based on the analogy to gravity, we 
may also derive the pseudo energy-momentum tensor of 
the <j> field with respect to the above metric g^ v 
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(d^){d v (l>)--g^ (5 P 0)(9^).(5) 



The associated energy density Tq of the light field 



n^ = - {[d^Y + (4 - v 2 ) [d x 



(6) 



contains negative parts beyond the horizon v 2 > A. Of 
course, this is precisely the reason why effects like Hawk- 
ing radiation are possible [11] . 

However, an energy density which is not bounded from 
below seems unphysical and typically indicates instabil- 
ities (already on the classical level). In order to avoid 
this problem, we may add an extra term which does not 
modify the linearized low-energy behavior of our model 



C^^C^a 2 (cl-v 2 f[d x ^ 
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but generates a positive definite energy density 



a(cl-v 2 )[d x <t>} 2 + ± 



(7) 



(8) 



In the exterior region c 2 > v 2 , the classical ground state 
is still given by (p = 0, but beyond the horizon A < v 2 , 
we have 2a(d x (p) = (v 2 — c 2 ) -1 / 2 . Thus, the classical 
ground state profile would not be differentiable at the 
horizon, i.e., the term [d^] 2 in the energy density, for 
example, would be ill-defined. This problem can be cured 
by adding another term (which again does not modify the 
low-energy behavior) and we finally arrive at the total 
Lagrangian of our toy model 



-full 



CI 



(9) 



The last term smoothens the classical ground state profile 
at the horizon and induces a super-luminal dispersion 
relation (uj + vk) 2 = A^k 2 + 2/3 2 fc 4 at large wavenumbers. 

Back-reaction The equation of motion of the light 
field can be derived from the Lagrangian above 

(d t + vd x )(d t + d x v)<f> = c 2 d 2 x cb + o(d%) , (io) 

where 0{d x ) denote the higher-order a and /3 terms we 
added for stability and regularity reasons. Similarly, the 
heavy field evolves according to 

i> - c\d 2 x ^ = V'W) ~ g[d t cj> + g^d x cj>]d x <t> + 0{dt) .(11) 

From the full set of equations, we see that the kink profile 
in Eq. @ together with </> = exactly solves the clas- 
sical equations of motion (though it is not the ground 



state). However, the impact of quantum fluctuations 
changes this picture: For 2irgipo > c<p, the kink acts as a 
black hole horizon and thus emits Hawking radiation. Of 
course, the energy /momentum given off must come from 
somewhere and hence this quantum effects should have 
some impact on the classical kink background. 

In order to estimate the quantum back-reaction, we 
quantize the fields <p — > cf> as well as ip — > tp and employ a 
mean-field expansion ip = ip c \ + Sip where ip c \ denotes the 
classical kink profile in Eq. @ and Sip as well as <p are 
supposed to be small (i.e., 4>,5ip <C ^ci)- Taking the ex- 
pectation value of Eq. (TTTj) and comparing it with Eq. (JSJ , 
we find that the lowest-order contributions of the quan- 
tum back-reaction force are just given by the expectation 
value of the pseudo energy-momentum tensor [12| 

[d 2 - c%d 2 x - V"^)] (SiP) » -g(f?) . (12) 
Remembering the covariant energy-momentum balance 

= {V~gTZ) - ~ T a " d v9aP = o , (13) 
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we find that (Tf) denotes the momentum density 
which varies with position in general. In contrast, the 
energy flux (Tq 1 ) measured with respect to the stationary 
frame is constant d x (T ( ]) = for a kink at rest. 

Fortunately, the expectation value (Tjf) can be calcu- 
lated analytically for a scalar field in 1+1 dimensions. 
In the Unruh state (which is the appropriate state for 
describing black- hole evaporation), one obtains [l3[ 



_ A-Vc^k 2 - [v'} 2 - jvv") - k 2 (c0 + v) 2 



<T°> = 



, (14) 



with 7 = 1 — AjA. and the effective surface gravity k 
determining the Hawking temperature 



1 { dv 



Hawking — 



2-7T 2tt \ dx 



(15) 



Note that (T®) calculated in the Unruh state is regular 
across black- hole horizon v = — c^, but singular at the 
white hole horizon v = +c^. (The Israel- Hartle-Hawking 
state would be regular at both horizons.) Far away from 
the kink/horizon v — > 0, we just get the usual thermal 
flux (Ti°) = -K 2 /(487rc ). 

The corrections induced by the quantum back-reaction 
can be visualized by incorporating them into an effective 
potential V e s via 



V^(iP) = V'(iP cl )-g(f?) 



(16) 



For the classical potential V(ip), all minima tp € 27r-0oZ 
occur at the same energy V = 0. However, the effec- 
tive potential V c r is distorted such that the central min- 
imum is lower than the next one describing the black 
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hole interior V c s(ijj — 0) < V c s(— 2tti/jq). In this sense, 
the exterior region is effectively energetically favorable 
and thus the horizon starts to move inwards, i.e., the 
black hole shrinks. Alternatively, the same result can 
be derived directly from Eq. (|12[) via classical time- 
dependent perturbation theory around the kink solution. 
The differential operator on the left-hand side of Eq. (fT2|) 
possesses a continuum of gapped propagating (delocal- 
ized) modes with to 2 > and one localized zero-mode 
oc 1/ cosh(^[x— Xkink]) with lo — 0, which just corresponds 
to a translation of the kink position [l4| ■ After expand- 
ing the source term —g(Ti) in Eq. (fT2f into these modes, 
the perturbations in the continuous spectrum oj 2 > just 
propagate away from the kink - whereas the spatial over- 
lap between — ^(Tj 3 ) and the zero-mode determines the 
acceleration ikink < of the kink position. 

Energy and Momentum In contrast to the fluid ana- 
logues for black holes (with a steady in- and out-flow 
of energy and momentum), for example, the kink con- 
sidered here represents a well localized object, which 
allows us to ask the question of where the force push- 
ing back the horizon comes from. In general, the con- 
tribution of the (f> field to the total energy-momentum 
tensor d^T^ = is different from the pseudo energy- 
momentum tensor V M T Ml/ = defined with respect to 
the effective metric (@|, which complicates the analysis 
[lit . Fortunately, these difficulties are absent in our toy 
model where the mixed components of both tensors coin- 
cide TJf = Ttf, The energy density T ° is given by Eq. © 
and the classical expression for the momentum flux den- 
sity just reads Tf = — Tq due to conformal invariance of 
the scalar field in 1+1 dimensions. Note, however, that 
the quantum expectation values differ due to the trace 
anomaly [13J. The energy flux density Tq 1 = <j)dC/d(f)' is 



given by Tq — cb[v(f) + (v 2 — A )<j>'\ and differs from the 
momentum density T® in Eqs. (fTTj) and (fT2|) for v ^ 0. 

Far away from the kink, we may estimate the above 
quantities by employing the geometric-optics approxima- 
tion and replacing c^> — > O and (j)' — > k. For solutions of 



the dispersion relation (SI + vky 



„2 1.2 
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responding to the outgoing Hawking radiation and its 
infalling partner particles, the energy density per normal- 
ized amplitude Tq = cSl 2 / (c — \ v\) changes its sign at the 
horizon, cf. Eq. ([6|). The energy flux density Tq = cSl 2 
is constant and positive everywhere (which is even true 
beyond the geometric-optics approximation). Note that 
SI is conserved as we are considering a quasi-stationary 
scenario. Thus, the total energy budget is balanced since 
the outgoing Hawking radiation carries away positive en- 
ergy, but the infalling partners have a negative energy. 

The momentum density T-j 3 = — cS! 2 /(c— \v\) 2 , on the 
other hand, turns out to be negative everywhere - or 
more precisely, far away from the kink, cf. the exact ex- 
pression (I14p with SI ~ k. Thus the momentum flux den- 
sity Tl = —cQ 2 /(c— \v\), i.e., the pressure, also changes 
sign at the horizon. (The trace anomaly vanishes in 



the asymptotic region v' = v" = far away from the 
kink where the geometric-optics approximation applies 
Tl = — Tq.) Consequently, while the Hawking particles 
carry away positive momentum and push back the kink, 
their infalling partner particles act in the opposite way 
and pull on the kink. In summary, the momentum is 
not balanced and thus the kink starts to move, i.e., the 
black-hole interior region shrinks. 

Thermodynamics The application of thermodynamic 
concepts to our toy model (in analogy to real black holes) 
presents some difficulties and ambiguities: Considering 
the heat capacity C = dE/dT, for example, we would 
associate T with the Hawking temperature (fTS)) . The 
variation of the internal energy dE 7 however, could be 
identified with the heat given off by the Hawking radia- 
tion dE = 5Q oc K 2 dt or with the change of the kinetic 
energy of the kink E = M cS xl ink /2 (for i^ ink < c|). 
Since the kink does not possess a conserved ADM mass, 
these quantities will be different in general. Either way, 
the heat capacity C = dE/dT could be positive as well 
as negative since the Hawking temperature can be in- 
creased dT > or decreased dT < by the quantum 
back-reaction of the evaporation process. There are sev- 
eral different effects: Due to the distortion of the effec- 
tive potential 14ff, the shape of the kink deviates from 
the classical profile ([2j) - This deviation is governed by 
the aforementioned continuum modes u> 2 > 0. Further- 
more, the kink starts to move - which is described by 
the zero-mode. The motion of the kink, in turn, implies 
a Doppler shift of the Hawking radiation. Finally, even 
in the rest frame of the kink, the position of the hori- 
zon Xh changes since the kink velocity cfckink effectively 
reduces the local frame-dragging speed v and therefore 
the surface gravity k = v'(xh) may change. As a result 
of all these effects, the heat capacity depends on many 
parameters (c^, c^, and gipo etc.) and may assume neg- 
ative as well as positive values. In order to demonstrate 
this sign ambiguity, let us consider the case G0 ^> cp for 
simplicity. In this limit, the continuum modes u) 2 > 
are very fast and hence the change of the shape of the 
kink can be neglected, i.e., the quantum back- reaction 
induces a rigid motion of the kink only. As another sim- 
plification, the transformation of the ?/;-field into the rest 
frame of the kink is just a Galilei transformation due to 
CtP 3> c^. The new horizon position is then simply deter- 
mined by v{xh) = — + ikink- Linearizing this equality 
together with k = v'(xh), we find that the variation 5k 
of the surface gravity induced by the acceleration of the 
kink <5±kink is determined by 5k — v" (xh)5%kink/ K- Since 
v"(xh) can be positive or negative (depending on the re- 
lation between and g4>o), the temperature measured 
in the rest frame of the kink could change in both di- 
rections. The temperature in the laboratory frame ac- 
quires an additional Doppler shift, which is given by 
5k = — K<5ikink/c</,. The relative strength of the two com- 
peting effects (Doppler shift and horizon displacement) 
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is given by c<pv" / (v') 2 , which can be above or below one. 
Ergo, both temperatures (in the kink frame and in the 
laboratory frame) may increase or decrease due to the 
back-reaction of Hawking radiation, i.e., the heat capac- 
ity can be positive or negative (or even infinite - at the 
turning point where ST = 0). 

Similar ambiguities apply to the entropy dS — dE/T. 
Choosing dE — SQ oc n 2 dt just reproduces the entropy 
balance of the Hawking radiation in the exterior region - 
which is of course indeed thermal. Inserting the kinetic 
energy E — AI c ffi 2 ink /2, on the other hand, we could 
violate the 2 nd law since the kink can be slowed down by 
incident coherent radiation (carrying zero entropy). 

Conclusions Modeling the black hole (horizon) by a 
stable topological defect in the form of a kink, we were 
able to derive the quantum back-reaction of the resulting 
evaporation process. It turns out that the kink/horizon is 
also pushed inwards as in a real black hole but, in contrast 
to the gravitational case, this back-reaction force is not 
caused by energy conservation but by momentum bal- 
ance. Energetically, the expansion of the horizon would 
be favorable because the minimum energy density in ex- 
terior region </> = ?/> = is far above l/(4a) 2 > the 
ground state in the interior region. Hence, going beyond 
the linear analysis performed here, one might suspect 
that the <fr field approaches its ground state via non-linear 
(quantum) instabilities until the evaporation stops. 

Further thermodynamical concepts such as heat ca- 
pacity or entropy (variation) cannot be defined unam- 
biguously and can have both signs - depending on the 
considered parameters [16(. Together with the results in 
0, our calculations and the energy-momentum consid- 
erations above suggest that Hawking radiation and the 
resulting back-reaction force "pushing" the horizon in- 
wards may be universal - whereas the heat capacity and 
the entropy concept strongly depend on the underlying 
structure (e.g., Einstein equations). Note that in the 
Israel-Hartle-Hawking state with the expectation value 
being (T°) = v(k 2 — [v'] 2 — ^vv") / {12ttc 2 ^ 2 ), the horizon 
is still pushed inwards - i.e., it does not correspond to 
the thermal equilibrium state for the combined system 
[kink in Eq. © plus <p field]. 
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